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In this paper, we present some asymptotic properties of the normalized inverse- 
^ ■ Gaussian process. In particular, when the concentration parameter is large, we estab- 

lish an analogue of the empirical functional central limit theorem, the strong law of 
large numbers and the Glivenko-Cantelli theorem for the normalized inverse-Gaussian 
process and its corresponding quantile process. We also derive a finite sum-representation 
that converges almost surely to the Ferguson and Klass representation of the normal- 
ized inverse-Gaussian process. This almost sure approximation can be used to simulate 
efficiently the normalized inverse-Gaussian process. 
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1 Introduction 

The objective of Bayesian nonparametric inference is to place a prior on the space of prob- 
ability measures. The Dirichlet process, formally introduced in Ferguson (1973), is consid- 
ered the first celebrated example on this space. Several alternatives of the Dirichlet process 
have been proposed in the literature. In this paper, we focus on one such prior, namely the 
normalized inverse-Gaussian process introduced by Lijoi, Mena and Priinster (2005). 

We begin by recalling the definition of the normalized inverse-Gaussian distribution. 
The random vector (Zi, . . . , Z m ) is said to have the normalized inverse-Gaussian distri- 
bution with parameters (71, ... , 7 m ), where 7i > for all i, if it has the joint probability 
density function 




in 



3/2 



X I s (zi, ...,z m ) 



(1.1) 



i=l 



where K is the modified Bessel function of the third type and § = {(21, . . . , z m ) : Z{ > 0, 
Y^iLi z i — l}- F° r more details about the modified bessel functions consult Abramowitz 



and Stegun (1972, Chapter 9). 
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Consider a space X with a o— algebra A of subsets of X. Let H be a fixed probability 
measure on (X, A) and a be a positive number. Following Lijoi, Mena and Priinster (2005), 
a random probability measure Pn, a = {Pn,a (A) } AeA is called a normalized inverse-Gaussian 
process on (X, A) with parameters a and H, if for any finite measurable partition Ax, ... , 
A m of X, the joint distribution of the vector (Ph,o,(Ai), . . . Pff )tt (A m )) has the normal- 
ized inverse-Gaussian distribution with parameter (aif(Ax), . . . aH(A m )). We assume 
that if H(Ai) = 0, then Pff, a (Aj) = with probability one. The normalized inverse- 
Gaussian process with parameters a and H is denoted by N-IGP(a, H), and we write 
P H , a ~X-IGP(a,H). 

One of the basic properties of the normalized inverse-Gaussian process is that for any 

A e A, 



E(P H>a (A)) = H(A) and Var(P H , a (A)) = H M\ (1.2) 

£(a) 



where here and throughout this paper 



«< a > = <13) 

andr(-2,#) = J^t^e-'dt. 

Furthermore, for any two disjoint sets Aj and Aj e .A, 

^F^A^P^A,-)) = ff(4)ff(A ^y (1.4) 

s \ a ) 

Observe that, for large a, £(a) ~ a (Abramowitz and Stegun, 1972, Formula 6.5.32, 
page 263), where we use the notation f(a) ps (7(a) if lim^^ f(a)/g(a) = 1. It follows 
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from (11.2b that H plays the role of the center of the process, while a can be viewed as the 
concentration parameter. The larger a is, the more likely it is that the realization of Pn, a is 
close to H . Specifically, for any fixed set A G A and e > 0, we have P H ^ a (A) A H(A) as 
a — > oo since 

ft{M . w>1 ,<«» as, 

Similar to the Dirichlet process, a series representation of the normalize inverse-Gaussian 
process can be easily derived from the Ferguson and Klass representation (1972). Specif- 
ically, let (Ei)i>x be a sequence of independent and identically distributed (i.i.d.) random 
variables with an exponential distribution with mean of 1 . Define 

Y i = E l + --- + E i . (1.6) 

Let (#i)i>i be a sequence of i.i.d. random variables with values in X and common dis- 
tribution H, independent of (r»)i>i. Then the normalize inverse-Gaussian process with 
parameter a and H can be expressed as a normalized series representation 

°° L -1 (T ) 

PhA-) = E ^L-ifr/ ^)' (L7) 



where 



POO 

L ( x ) = " / e -*/ 2 t- 3 / 2 tZt, for a: > 0, 



(1.8) 



and 6x denotes the Dirac measure at X (i.e. 8x(B) = 1 if X E B and otherwise). Note 
that, working with (11.7b is difficult in practice because no closed form for the inverse of the 
Levy measure (11.8b exists. Moreover, to determine the random weights in (11.7b an infinite 
sum must be computed. 
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This paper is organized as follows. In Section 2, we study the weak convergence of the 
centered and scaled process 

as a — > oo. Note that, since £(a) ~ a, it is possible to replace the normalizing coefficient 
£(a) by a. Therefore, for simplicity, we focus on the process 

D Hta (-) = ^(P HA (.)-H(-)). (1.9) 

In Section 3, we derive the limiting process for the quantile process 

QhA-) = ^(p^(-)-h-H-)), (i.io) 

as a — > oo, where, in general, the inverse of a distribution function F is defined by 

F~ l (t) = inf {x : F(x) >t}, < t < 1. 

The strong law of large numbers and the Glivenko-Cantelli theorem for the normalized 
inverse-Gaussian process are discussed in Section 4. In Section 5, we derive a finite sum- 
representation which converges almost surely to the Ferguson and Klass representation of 
the normalized inverse-Gaussian process. This new representation provides a simple, yet 
efficient procedure, to sample the the normalized inverse-Gaussian process. 
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2 Asymptotic Properties of the Normalized Inverse-Gaussian 
Process 

In this section, we study the weak convergence of the process D H a defined in (11.91) for 
large values of a. Let y be a collection of Borel sets in R and H be a probability measure 
on R. We recall the definition of a Brownian bridge indexed by ,5? . A Gaussian process 
{Bh(S) : S G J^} is called a Brownian bridge with parameter measure H if E [Bh{S)] = 
for any S G 5? and 

Cov (Bji(Si), B H (Sj)) = H(Si n Sj) - H(Si)H(Sj) (2.1) 

for any S>, Sj G S? (Kim and Bickel, 2003). 

The next lemma gives the limiting distribution of the process (11.91 ) for any finite Borel 
set Si, ... , S m G y, as a — > oo. The proof of the lemma for m = 2 is given in the 
appendix and can be generalized easily to the case of arbitrary m. 

Lemma 1. Let D H a be defined by (fi.9D . For any fixed sets Si, . . . , S m in 5? we have 

(D Hia (Si),D Hia (S 2 ), D Hta (S m )) A (B H (Si), B H (S 2 ), B H (S m )) , 

as a — >■ oo, where B H is the Brownian bridge with parameter H. 

The following theorems show that the process Du, a defined by (11.91) converges to the 
process B H on D[— oo,oo] with respect to the Skorokhod topology, where D\— oo, oo] 
is the space of cadlag functions (right continuous with left limits) on [—00,00]. Right 
continuity at —00 can be achieved by setting D H>a (—oo) = D Hta (—oo) = 0; the left limit 
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at oo also equals zero, the natural value of DhA°°) and D Ha {oo). For more details, 
consult Pollard (1984, Chapter 5). If X and (X a ) a>0 are random variables with values in a 
metric space M, we say that (X a ) a converges in distribution to X as a — > oo (and we write 
X a -4 X) if for any sequence (a n ) n converging to oo, X an converges in distribution to X. 

Theorem 1. We have, as a — » oo, 

DhA-) = (PhA-) ~ #(•)) ^ b h(-) (2.2) 

on D[— oo, oo] with respect to Skorokhod topology, where B H is the Brownian bridge with 
parameter measure H. 

Remark 1. In the proof of Lemma Q3 the finite-dimensional convergence is in fact con- 
vergence in total variation, which is stronger than convergence in distribution (Billingsley 
1999, page 29). 

Remark 2. Lo (1987) obtained a result similar to that given in Theorem Q] for the Dirichlet 
process to establish asymptotic validity of the Bayesian bootstrap. An interesting general- 
ization of Lo (1987) to the two-parameter Poisson-Dirichlet process was obtained by James 
(2008). In both papers, proofs are based on constructing the distributional identities (Propo- 
sition 4.1, James, 2008)). Establishing analogous distributional identity for the normalized 
inverse-Gaussian process does not seem to be trivial. 

Remark 3. Sethuraman and Tiwari (1982) studied the convergence and tightness of the 
Dirichlet process as the parameters are allowed to converge in a certain sense. Analogous 
results for the two-parameter Poisson-Dirichlet process (Pitman and Yor, 1997; Ishwaran 
and James, 2001) and the normalized inverse-Gaussian process (PH, a ) a>0 follows straight- 
forwardly by applying the technique used in the proof of the tightness part (i.e. condition 
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d73b ) of Theorem Q] 



3 Asymptotic Properties of the Normalized Inverse-Gaussian 
Quantile Process 

Similar to the frequentist asymptotic theory, in this section we establish large sample theory 
for the normalized inverse-Gaussian quantile process. 

Corollary 1. Let < p < q < 1, and H bea continuous function with positive derivative h 
on the interval [H~ 1 (p) — e, H~ l (q) + e] for some e > 0. Let A be the Lebesgue measure on 
[0, 1]. Let Qh,q be the normalized inverse-Gaussian process defined in (fl.iOj) . As a — > oo, 
we have 

in D\p,q\. That is, the limiting process is a Gaussian process with zero-mean and covari- 
ance function 

^ x X(s At) - X(s)X(t) 

Proof. We only proof part (i). Part (ii) follows similarly. By Theorem \T\ (i) the pro- 
cess \fa(Pu,a — H) converges in distribution to the process Bh = B\{H) = B\ o H. 
Notice that almost all sample paths of the limiting process are continuous on the inter- 
val [H~ l {p] — e, H~ 1 (q) + e] . By Lemma 3.9.23 page 386 of van der Vaart and Wellner 
(1996), the inverse map H i— )■ H" 1 is Hadamard tangentially differentiable at H to the 
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subspace of functions that are continuous on this interval. By the functional delta method 
(Theorem 3.9.4 page 374 of van der Vaart and Wellner (1996)) we have 



Remark 4. A similar result to Corollary \T\ for the two-parameter Poisson-Dirichlet process 
can be obtain by applying Theorem 4.1 and Theorem 4.2 of James (2008). 
Remark 5. Similar to Remark 1 of Bickel and Freedman (1981), if H~ 1 (0+) > — oo and 
iJ _1 (l) < oo and h is continuous on [H~ 1 (0+], the conclusion of the corollary 

holds in D [H~ 1 (0+), iJ _1 (l)] . For example, if H is a uniform distribution on [0, 1], then 
the convergence holds in D[0, 1]. 

The next example is a direct application of Corollary [Q 
Example 1. In this example we derive the asymptotic distribution for the median and the 
interquantile range for the the normalized inverse-Gaussian process. Let Q X H a , Q 2 H a and 
Q Z H a be the first, the second (median) and the third quartiles of Ph,q (i- e - ^#1(0-25) = 

Qh,* ^K - 5 ) = Q 2 H,a and P hU - 75 ) = Qh,o)- Let ?i. 12 and q 3 be the first, the second 
(median) and the third quartiles of H. From CorollaryCQ after some simple calculations, the 
asymptotic distribution of the median and the interquantile range are given, respectively, 



QhA-) ^ 



B x oHoH^(- 



h{H-*{.)) 



h{H-\-)) 



in D[p, q}. This completes the proof of the corollary. 



□ 



by: 
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and 

V-aWn - fa, - A n (o, ^ + ^ - jj^) , 

where h = H' and JQ-R = — <5}/ a - Note that, the asymptotic distributions of the 
median and the interquantile range for the normalized inverse-Gaussian process coincide 
with that of the sample median and the sample interquartile range (DasGupta, 2008, page 
93). 

4 Glivenko-Cantelli Theorem for the Normalized Inverse- 
Gaussian Process 

In this section, we show that an analogue of the empirical strong law of large numbers 
and the empirical Glivenko-Cantelli theorem continue to hold for the normalized inverse- 
Gaussian process. 

Theorem 2. Let P H a ~ N-IGP(a, H). Assume that a = n 2 c, for a fixed positive number 
c. Then as n — > oo, 

Ph,ti 2 c(A) a A H(A), 

for any measurable subset A of X. 
Proof. For any e > 0, by (14.11) . we have 

P r {M A)-H { A)l> e) < H ^-f A ^\ (4.1) 
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where £(n 2 c) is defined by (11.31) . Note that 

n-*oo 1/n C 

(Abramowitz and Stegun, 1972, Formula 6.5.32, page 263). Since the series Yl^Li V^ 2 
converges, it follows by the Limit Comparison Test that the series Yl^=i £ (^ 2 c) is a ls° 
convergent. Thus, 

oo 

^Pr{|P Hi „ 2c (A) - > e} < oo. 

n=l 

Therefore, by the first Borel-Cantelli Lemma, the proof follows. □ 

The proof of the next theorem follows by arguments similar to that given in the proof 
of the Glivenko-Cantelli theorem for the empirical process. See, for example, Billingsely 
(1995, Theorem 20.6). 

Theorem 3. Let Pu,a ~ N-IGP(a, H). Assume that a = n 2 c, for a fixed positive number 
c. Then 

sup\P H , n 2 c (x)-H(x)\ a -$0, (4.2) 

as n — y oo. 

Remark 6. Similar to the normalized inverse-Gaussian process, a strong law of large 
numbers and a Glivenko-Cantelli theorem can also be established for the two-parameter 
Poisson-Dirichlet process. 
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5 Monotonically Decreasing Approximation to the Inverse 
Gaussian Process 

In this section, we derive a finite sum representation which converges almost surely to the 
Ferguson and Klass sum representation of the normalized inverse-Gaussian process. We 
mimic the approach developed recently by Zarepour and Al Labadi (2012) for the Dirichlet 
process. Let X n be a random variable with distribution inverse-Gaussian with parameter 
a/n and 1 (see equation (3) of Lijoi, Mena and Priinster (2005) for the density of the 
inverse-Gaussian distribution). Define 

G n (x) = Pr(X n > x) = -^—r 3 / 2 exp{-U^-+t)+-)dt. (5.1) 

The following proposition describes properties of G n (x) that will be used later in the 
paper. 

Proposition 1. Forx > 0, the function G n (x) defined in (f5.il) has the following properties 

as n — > oo: 

(i) nG n (x) L(x), 

(ii) G-\l) -> L~\x), 
where L is defined in (2.3). 

Proof. To prove (i), since e a ^ n < e a and exp j — | {^- t + ^) } — 1> me integrand in 
nG n (x) is bounded by ae a t~ 3 ^ 2 j\p2m, which is integrable for any x > 0. Hence, the 
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dominated convergence theorem applies and we have 



a 



POD 

/ e-'lH-*' 2 dt = L{x) 



nG n (x) = 



dt 



a 



V2tt 



t~ 



V2tt 



J x 



To prove (ii), notice that the left hand side of (i) is a sequence of monotone functions 
converging to a continuous monotone function for every x > (Haan-de and Ferreira 2006, 



PropositionCQgives a simple procedure for an approximate evaluation of both L(x) and 
L~ l (x) for any x > 0. For computational simplicity, a more convenient approximation is 
presented in the following Corollary. The proof follows straightforwardly by taking x = Tj 
in Proposition 1 and the fact that we have T n+1 /n ^4' 1 as n — > oo (strong law of large 
numbers). 

Corollary 2. For a fixed i, as n — > oo, we have: 



Remark 7. The utility of Corollary [2] stems from the fact that all values of G n 1 (r i /r n+ i) 
are nonzero for i < n. This is not the case when working with G~ 1 (r i /n). 

The following lemma provides a finite sum representation which converges, almost 
surely, to the Ferguson and Klass (1972) sum-representation for the inverse-Gaussian pro- 
cess. The proof of the lemma is similar to that of Lemma 2 in Zarepour and Al Labadi 
(2012). Hence, it is omitted. 

Lemma 2. If (0;)i>i is a sequence of i.i.d. random variables with common distribution H, 



page 5). Thus, (i) is equivalent to G r n 1 (x/n) — >■ N 



□ 
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independent of (Vi)i>i, then as n — > oo 

i=l ^ n+l / 

Here, r i; iV(x), and G n (x), are defined in (U.6D . (fi.8D . and 0. J D . respectively. 

By normalizing the finite sum in (15.21) , it is possible to obtain a sum representation that 
converges almost surely to Ferguson and Klass representation of the normalized inverse- 
Gaussian process. This important result is stated formally in the next theorem. 

Theorem 4. Let {9i)i>\ he a sequence of i.i.d. random variables with values in X and 
common distribution H, independent of (Ti)i>i, then as n — > 00 



( r< \ 


\Tn- 


-1 J 


w 






r n +i 



(5.3) 



Here Ti, L(x), and G n {x), are defined in (fj.6D . (U.8D . and (f5.JD . respectively. 

Remark 7. For any 1 < i < n, rj/r n+1 < T i+ i/T n+ i almost surely. Since G~ l is a 
decreasing function, we have G^ 1 (Ti/T n+1 ) > G" 1 (T i+1 /T n+1 ) almost surely. That is, 
the weights of the new representation given in Theorem 0] decrease monotonically for any 
fixed positive integer n. As demonstrated in Zarepour and Al Labadi (2012) in the case of 
the Dirichlet process, we anticipate that this new representation will yield highly accurate 
approximations to the normalized inverse-Gaussian process. 
Remark 7. For a of Theorem |4] we can write 



71 

d 



i=l 



Al Labadi and Zarepour: Asymptotic Properties of the Normalized Inverse-Gaussian Process\5 

where pi )Tl , . . . ,p n , n ~ N-IG (a/n, . . . , a/n), = means have the same distribution and the 
N-IG distribution is given by (11.11) . Therefore a similar result to Theorem 2 of Ishwaran 
and Zarepour for the normalized inverse-Gaussian process follows immediately. 

6 Concluding Remarks 

The approach used in this paper can be applied to similar processes with tractable finite 
dimensional distributions. On the other hand, when the finite dimensional distribution is 
unknown, one may follow the approach of James (2008). However, applying this approach 
requires constructing distributional identities, which may be difficult in some cases. 

One can use the results obtained in this paper to derive asymptotic properties of any 
Hadamard-differentiable functional of the N-IGP(a, H) as a — > oo. For different applica- 
tions in statistics we refer the reader to van der Vaart and Wellner (1996, Section 3.9) and 
Lo (1987). Moreover, it is possible to extend the results found in this paper to the case 
when the base measure H is a multivariate cumulative distribution function. The result of 
Bickle and Wichura (1972) can be employed in the proof. 
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Appendix 

Proof of Lemma 1 for m = 2: Let Si and S 2 be any two intervals in K. Without loss 
of generality, we assume that Si D £2 = 0- The general case when Si and S 2 are not 
necessarily disjoint follows from the continuous mapping theorem. 



where the N-IG distribution is given by (11.11) . For notational simplicity, set Xi = PH, a (Si), 
U = H(Si) and D t = y/a(Xi — l { ) for i = 1,2. Thus, the joint density function of Xi and 
X 2 is: 



Note that 



(PH,a(Sl), PH,a(S 2 ), 1 



N-lG(aH(Si),aH(S 2 ), 
a(l-H(Si)-H(S 2 ))), 



fxi,x 2 ( x i, x v 



e a a%l 2 (l -h -l 2 ) 

2 1 /%3/2 



X X 




(1 - xi - x 2 )-^ 




a 



1 / 2 e a / 1 / 2 (l -h- l 2 ) 



x x 




{i-xi-x 2 y^ 



21/2^3/2 




The joint probability density function of Di = ^fa [X x — li) and D 2 = y/a (X 2 — l 2 ) 
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is: 



a 1 / 2 e a ZiZ 2 (l -h- l 2 ) 

2V27T3/2 X 

x (yx/v^ + Z!)- 372 (y 2 /yfi + h)~ 3/2 
(l - Ui/Va- h ~ Vilsfa- l 2 ) 3/2 

(. « 



xK. 



-3/2 



72 
'2 



\ 1/2n 

(l-/l-/ 2 ) 2 



1 - Z/i/v 7 ^ - h - Vi/y/a - i 2 



72 
'1 



72 
'2 



Vi/y/a + h V2/Va + k 

(l-/i-/ 2 ) 2 
1 - Vi/y/a - h - U2/y/a - l 2 



-3/4 



By Scheffe's theorem (Billingsely 1999, page 29), it is enough to show that: 



fD 1 ,D 2 {Vl,V2) -> /(yi,y2) 



2tt|S| 1 /2 



exp{-(^ 2 )£-% l2 / 2 ) T /2}, (7.1) 



where £ 



Zi(l-Zi) -/i/ 2 
-W 2 z 2 (1 - i 2 ) 



Since, for large z and fixed z/, K v {z) ~ ^/tt/2z 1 ^ 2 e 2 (Abramowitz and Stegun, 1972, 
Formula 9.7.2, page 378), where we use the notation /(V) « gf(^) if lirn^oo f(z)/g(z) = 
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1, we get: 



lim f Dl ,D 2 (yi,y 2 ) 



Notice that, 



lim 

a— >ao 



hhiX-k-h) 



2tt 



x 



r 2 

'"2 



(l-h-k 



yyi/Va + k y 2 /\/a + l 2 1 - yi/y/a- k - y 2 /Va- k 

x (yi/V^ + hT 3/2 fa/VZ+h)'*'* 
(l-yi/y/a-l 1 -y 2 /y/a-l 2 ) 3/2 



x exp a 



72 
'2 



\yi/\fa + l\ y2/Va + h 

1/2, 



(1-h- k) 



1 - yi/Va - h - Vz/Va - h 



hl 2 {l-h-k 



x 



'2 



(i-h-ky 



-i 



\yi/Va + h y2/Va + k i-yi/Va-h-y2/Va-l 2 
x ( yi /VE + h)~ 3/2 (y 2 /V^ + ky V2 (1 - yi/V^-h ~ y*l\fc- h)~ 3/2 



converges to 1/ (2n y/a n a 2 2(l — P12)) > where 



(Tn = h(l - h), a 22 = l 2 {l-l 2 ), P12 



(l-h)(l-l 2 ) 



To prove the lemma, it remains to show that 



a 1 



7 2 
'1 



72 
f 2 



(l-k-hY 



l/2> 



yi/Va + k Vz/y/a + h l-yilVa-li-y 2 /sfa-l 2 
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converges to 

The last argument follows straightforwardly form the L'Hospital's rule. □ 

Proof of Theorem 1: Let (a n ) be an arbitrary sequence such that a n — > oo. To simplify 
the notation, in the argument below, we omit writing the index n of a n . Assume first that 
H(t) = X(t) = t (i.e. A is the Lebesgue measure on [0, 1]). Thus the process (12.21 ) reduces 
to 

D x ^{t) = ^{P Ka {t)-t). 

To prove the theorem, we use Lemma[T|and Theorem 13.5 of Billingsley (1999). Therefore, 
we only need to show that for any < ti < t < t 2 < 1, 

< \F(t 2 ) - F^f 1 , (7.2) 

for some (3 > 0, a > 1/2, and a nondecreasing continuous function F on [0, 1]. Take 
a = ft = 1 and F(t) = t. We show that 

E [(D A ,„(t) - D^t,)) 2 (D x , a (t 2 ) - D^ a {t)f] < (t 2 - t,) 2 . (7.3) 

Observe that 



1 

2(1 - p\ 2 ) 



E 



D x , a (t) - D x , a {ti)\ 2f> ~ Dx, a (t)\ 



2/3 



D x ,a{t) ~ Dx,a{h) = Da,b((<i,<]) and D x , a (t 2 ) - D Aj0 (t) = D x>a ({t,t 2 ]). 
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Thus, the expectation in the right hand side of ( 17.31) is equal to 

e(a)E - A((t 1? t])} 2 {P A , a ((M 2 ]) - A((M 2 ])} 2 ] , (7.4) 

where A((£, t 2 ]) = t 2 — t and A((ii, t\) —t — t\. Expanding the expression 

{Px, a {{hA) - K(ti,t])} 2 {Px, a ((tM) - mM)} 2 (7-5) 

gives 

Pla((h,t])Pla((t, t 2 \) - 2A((t, t 2 ])Pl a (( tll t])P X , a ((t, t 2 ]) 

+ A 2 ((t,t 2 ])P A 2 a ((t 1 ,t])-2A((t 1 ,t])P A , a ((t 1) t])P A 2 a ((t,t 2 ]) 

+ 4A((ti, t])A((t, t 2 ])PA,a((tl, t])P A ,a((t, t 2 ]) - 2A((ti, t])A 2 ((t, t 2 ])P A ,a((tl, *]) 

+ A 2 ((t 1 ,t])P A 2 a ((M 2 ])-2A 2 ^ 

Using the fact that P A , a (-) is a probability measure and < t\ < t < t 2 < 1, the expression 
displayed in (17.51 ) is less than or equal to 

n 2 a ((^i,t])P A 2 a ((t,t 2 ]) + A 2 ((t,t 2 ])P A 2 a ((t 1 ,t]) 
+4A((t 1 ,t])A((t,t 2 ])P A , a ((t 1 ,t])P A , a ((t,t 2 ]) + A 2 ((t 1 ,t])P A 2 a ((t,t 2 ]) 
+A 2 ((t 1 ,t])A 2 ((M 2 ]) 

< P A , a ((ti, t])P A ,a((*, * 2 ]) + A((t, t 2 ])P A ,a((ti, *]) 

+4A((ti, t]) A((t, t 2 ]) + A((t 1; t])P A , a ((t, * 2 ]) + A((t, t 2 })X((t, t 2 }). 
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By (ITT2T ) and ( |TT4b we obtain 

E[{P Aia ((t 1 ,t])-A((t 1 ,t])} 2 {P A , a ((t ! t 2 ])-A((t,t 2 ])} 2 ] < 

a 

A((ti,t])A((t,t 2 ]). (7.6) 

Thus, using (17.41) and (17.61) . we have 

£ [(£>*,„(*) " ^(tl)) 2 pA,a(* 2 ) - D^ a {t)f] = 



< 

for < t\ < t < t2 < 1. This proves the theorem in the case when = £, i.e. H is 
the uniform distribution. Observe that, the quantile function _ff~ 1 (s) = inf {t : H(t) > s} 
has the property: H^ 1 (s) < t if and only if s < iP(t). If [/, is uniformly distributed over 
[0, 1], then H~ l (Ui) has distribution H. Thus, we can use the representation a, = 
where (C/f)i>i is a sequence of i.i.d. random variables with uniform distribution on [0, 1], 
to have: 

P H ,ait) = P Ka{H{t)) and D Hja {t) = D Ka {H{t)) = D Ka o Jf (*), t G R, 

where P AjQ is the normalized inverse-Gaussian process with concentration parameter a and 
Lebesgue base measure A on [0, 1]. From the uniform case, which was already treated, 
we have D A)0 (-) = y/E (P A , a (-) - A(-)) A P A (-)- Define : Z>[0, 1] P[-oo, oo] by 
(iffx)(t) = x(H(t)). Since the function ^ is uniformly continuous (Billingsley 1999, page 
150; Pollard, 1984, page 97), it follows, from the continuous mapping theorem and the fact 



8a- 
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a 3 




8a- 
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a 3 




8a- 


1 



A(ti,t]A(M 2 ] 
{t - h) (t 2 - t) 
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that D Xa A B x , that D H>a = $(D X;a ) -> ^(B x ) = B H . This completes the proof of the 
theorem. □ 



